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0.  Introduction 

Methods  of  estimating  the  scale  parameter  of  exponential  life 
distributions  are  generally  based  on  complete  information  of  the  failure 
times,  with  or  without  censoring.  More  specifically,  the  estimation  is 
based  on  life  testing  of  n  identical  systems  which  are  subjected  to 
similar  experimental  conditions.  The  failure  times  of  the  systems  are 
recorded,  and  the  experiment  is  censored  either  at  the  r-th  failure 
(1  <  r  <  n  )  ,  which  is  called  censoring  of  Type  II,  or  at  a  fixed  time, 

tQ  ,  which  is  censoring  of  type  I.  In  some  cases  the  minimum  of  tQ  and 

the  r-th  failure  is  used.  The  reader  is  referred  to  any  one  of  the  text¬ 
books  on  statistical  methods  in  reliability  for  the  details  of  the 
estimation  procedures  (see  for  example  W.  Nelson  [5]).  The  present 
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paper  has  been  motivated  by  two  papers  of  Epstein  et  al  [2,3],  which  discuss 

photodynamic  bioassays  for  estimating  the  potency  of  benzo-solubale  organic 

extracts  from  particulates  in  air  samples.  In  these  experiments  various 

dosages  of  organic  extracts  and  of  their  derived  fractions  (basic,  acidic, 

aliphatic,  aromatic,  oxygenated  and  water  soluble)  were  applied  on  n  ■  30 

living  organisms  (Paramicium  coudatum)  under  ultraviolet  irradiation. 

In  each  application  the  measured  biological  response  was  the  time  required 

(in  minutes)  to  immobilize  90%  of  the  cells.  This  response  variable  was 

called  the  lethal  time  90  (LT90) .  Each  trial  was  censored  at  t  *  90 

o 

minutes,  and  the  number  of  cells  still  living  at  the  time  of  censoring 
was  not  recorded.  Empirical  fitting  of  a  dosage-response  curve  was  done 
by  Epstein  et  al  [2].  Serious  questions  arise  concerning  the  validity 
of  such  an  empirical  analysis,  which  is  performed  without  sufficient 
theoretical  justification.  The  present  paper  is  devoted  to  a  theoretical 
analysis  of  two  possible  estimators  of  a  scale  parameter,  when  the  data 
are  the  censored  r-th  order  statistic.  More  specifically,  let  X 

n,r 

denote  the  r-th  order  statistic  of  a  sample  of  n  i.i.d.  random 

variables  from  an  exponential  distribution  with  mean  6  .  Let 

*  * 

X  =  min(X  ,  t  )  .  A  random  sample  of  m  observations  on  X 
n,r  n,r  o  nfr 

is  available.  In  Section  1  we  present  the  model  and  define  some  basic 

notions*  In  addition  we  provide  formulae  of  the  expected  value  and 
* 

variance  of  X  .On  the  basis  of  these  formulae  ve  study  in 
n,r  J 

Section  2  the  structure  and  properties  of  a  moment-equation  estimator 

r  *  i 

(MEE)  of  6  .  This  is  the  value  of  0  for  which  EQ{X  }  is  equal 

o  n,r 

* 

to  the  sample  mean  of  the  m  observations  on  X  .We  show  that 

n,r 

if  at  least  one  observation  is  not  censored  then  the  MEE  is  unique. 

We  further  show  that  as  m  grows,  the  asymptotic  distribution  of  the 
MEE  is  normal.  Formulae  for  the  asymptotic  bias  and  asymptotic  variance 
of  the  MEE  are  also  derived.  In  Section  3  we  study  the  structure  of 
the  maximum  likelihood  estimator  (MLE)  of  6  and  show  that  it  is  unique, 
provided  at  least  one  observation  is  uncensored.  The  Fisher  Information 
Function  and  the  asymptotic  normal  distribution  of  the  MLE  are  derived 


-  2  - 


T-495 


in  Section  4.  In  Section  S  ve  discuss  the  asymptotic  relative  efficiency 
of  the  MEE  compared  to  the  MLE,  and  provide  some  simulation  results  which 
illustrate  the  performance  of  these  estimators  in  samples  of  size 
m  *  20(10)50.  The  results  of  the  present  study  can  be  applied  not  only 
to  the  analysis  of  the  photodynamic  bioassays  data  of  Epstein  et  al  [2,3] 
but  also  to  reliability  studies,  in  which  the  observed  failure  times  of 
systems  are  based  on  censored  r-th  order  statistics.  This  is  the  case, 
for  example,  when  n  identical  components  within  a  system  function 
independently,  and  the  system  fails  when  the  r-th  component  failures 
happen.  The  times  of  the  component  failures  are  not  available. 

1.  The  Model  and  Moments  of  the  Censored  Observations 

Consider  a  sequence  of  independent  replicas  of  an  experiment  in 

which  n  systems  are  subjected  tQ  life  testing  (possibly  accelerated). 

Each  experiment  is  terminated  either  at  the  r-th  failure  (1  <  r  <  n)  , 

or  at  a  fixed  time  point,  tQ,  whichever  comes  first.  The  recorded  random 

variable  is  the  time-censored  r-th  order  statistic.  We  further  assume 

that  the  life  distribution  of  the  individual  systems  is  an  exponential 

distribution,  i.e.,  F(x/0)  *  1  -  exp  {-x/0}  ,  0  <  x  <  00  and  0  <  0  <  « 

Let  X  denote  the  r-th  order  statistic  in  a  sample  of  n  i.i.d. 
n,r 

X 

random  variables,  having  a  common  distribution  F(-r-)  .  The  observed 

®  * 

random  variable  in  the  time-censored  experiment  ic  X  *  min(t  ,  X  )  • 

o  n,r 

Since  the  life-time  distribution  F(^-)  is  exponential,  the  p.d.f.  of 

X  is  (see  David  [1]) 
n ,  r 

(1.1)  fn^r(x;0)  -  r(”)  |  exp{-(n-r+l)  (l-exp{-  |})r_1  , 

for  0  <  x  <  ®>  .  Let  U  *  X  /0  be  the  standardized  r-th  order 

n,r  n,r 

statistic,  and  let  f(u)  denote  the  standard  density,  i.e., 

fn  •  Let  ®  e»  O  designate  the  c.d.f .  of 

n,r  t?  n,r  o  n,r  o  ° 

* 

X  •  This  c.d.f.  is  absolutely  continuous  on  [0,t  )  and  has  a  lump 

o 

point  at  t  ,  i.e., 
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(1.2)  G  (x;  0,  t  ) 
n,  r  o 


,  x  <  0 


/  fn>r(x';  6)  dx' 


,  0  <  x  <  t 


*  x  >  fco  * 


The  height  of  the  jump  of  this  c.d.f.  at  x  «  t  is 


(1.3)  PA{X  >  t  )  »  B(r-1;  n,  F(-q))  , 

y  n, r  -  o  u 


where  B(j;  n,  p)  designates  the  c.d.f.  of  the  Binomial  distribution  with 


parameters  (n,  p)  .  We  derive  now  formulae  for  the  expected  value  and 
variance  of  the  censored  variable  X  .  It  is  well  known  (see  David  [1]) 


*  r_1  i 

if  t  *  00  ,  (no  censoring)  then  E  {X  |0,  t  *  »}  *  0  I  — t 

°  j*o 


*  2  1 
and  V{X  1 6,  t  *  »}  *  0  Z  - , 

°  j-o  (n-j)2 


.  This  is  due  to  the  fact  that 


the  r-th  order  statistic  in  the  exponential  case  can  be  represented  as  a 
sum  of  r  independent  random  variables,  having  exponential  distributions 
with  means  6/ (n-j)  ,  j-0, ...,r-l  .  In  the  time  censored  case  one  obtains 


Lemma  1. 1.  The  first  two  moments  of  X  in  the  time  censored  case  are: 


.  r-1  t 

(1.4)  E{X  |6,  t  }  -  0  I  0-B(j;  n,  F(^))] 

j*o  J 


(1.5)  E(X*2 1 6 ,  t  }  =  262  (V  -±r  Z  ^  [1  -  B(k;n,  F(^|))] 

(j*o  J  k=o 

t  r-1  t  \ 

-f  -3zmo  ^  B(j! »• F(-f»  | 
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r*T*' 


According  to  (1.2)  , 


(1.6)  E{X*|e,  t  }  -  /  x  f  (x;0)  dx  +  t  B(r-1;  n,  F(-|)) 

o  n,r  o  0 


*2  fP  o  o  t 

(1.7)  E{X  |e,  t  }  *  /  X  f  (x;0)  dx  +  t*  B(r-1;  n,  F(-£)) 
o  n  j  r  o  0 


Let  n  =  t  /0  then 


(1.8)  /  x  f  (x;0)  dx  *  0  /uf  (u)  du  =  0  yf^Ol)  . 

n | i  ^  u ) l  n,r 


Moreover,  by  interchanging  the  order  of  integration  we  obtain 


L.9)  u(1)(rtf  -  Jf  </dt)  f  (u)  du 
n.r  *  *  n,r 


n  n 


Furthermore, 


/  (/  fn,r(u)  du)  dt 

n  t  t 

/  [B(r-1;  n,  F(|))  -  B(r-1;  n,  F (-£))]  dt 


■  J  B(r-1;  n,  l-e_t)  dt  -  n  B(r-1;  n,  l-e~n) 


d t  -  T.  f  -(n-j)tf 


(1.10)  /B(r-1;  n,  l-e  fc)  dt  -  Z  (")  /  (l-e"V  e' 

o  j»o  J  o 
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j-o  'y  <j+i)  (^) ;  ' 


E  -4  u-»(j;  n,  l-e"n>]  . 

J-o  n'J 


Furthermore, 


Finally,  substituting  (1.9)  and  (1.10)  in  (1.6)  we  obtain  (1.4).  Notice 

that  lim  B  (j;  n,  i-e-ri)  *  0  for  each  j  <  n  .  Hence,  the  well  known 
n-K*> 

formula  of  E{X  ] 0 ,  t  =  “>}  is  a  limiting  case  of  (1.4).  In  a  similar 
fashion,  we  derive  the  formula 

(1.11)  P^Ol)  *  /  u2f  „(u)  du 
n.r  ^  n.r 

o 


2/</  tdt)  fQ  r(u)  du 


2  J  t  ( J  f^  r(u)du)dt 


2  J  t  B(r-1;  n,  l-e_t)dt  -  n2B(r-l;  n,  l-e”*1) 
o 


r-1  „  n 


(1.12)  /  t  B(r-1;  n,  l-e_t)dt  -  l  (”)  /  t  (l-e-t)j  e"t(n_;i)dt 

o  j*o  ^  o 

r-i  n 

-  Z  ~ -  /  t  f  .+1(t)dt 

j-o  n"j  o  ,J+1 

-  V  ^j[iio  ih  (»-»»;  o,  l-e'n))  -  n  »(Ji  »,  1-e-";] 


-  6  - 


A  *.  A  _%  A  *-  Jp 


%  ^  k> W  /V-VkVW  kV  » YV> 


r^r^r^rflr  wrt'rt‘,*r/r«'/r^,T<Jr7r .■»  t 


Substituting  (1.11)  and  (1.12)  in  (1.7)  with  the  proper  scale  factor,  we 
obtain  formula  (1.5). 

(Q.E.D.) 


From  formulae  (1.4)  and  (1.5)  one  obtains  immediately  the  following  formula 

* 

for  the  variance  of  X  ,  namely: 

n,r 


(1.13)  o2(0;  t  )  =  V{X*  ;  0,  t  } 

o  n,  r  o 


=  0 


(  r“l  t  J  ,  t 

2.z  ^  tz  S*  [1"B(k;  “■  f(t))] 

l  J=o  J  k=o 


-  *-t  y  ^  -•  f(-t»  - 

j=o  J 

2 

(X  »>  F(^f  >))  |  • 


2  2  1 

One  can  easily  check  that  lim  a  (0;  t  )  *  0  Z  - r 

tQ->co  °  j=o  (n-j) 

In  the  following  section  we  apply  the  above  formulae  for  the  study 
of  the  moment-equation  estimator,  MEE. 


2.  The  Moment  Equation  Estimator  and  its  Asymptotic  Properties 

In  the  case  of  no  time  censoring,  i.e.,  when  t  *  00  ,  an 

unbiased  estimator  of  6  is  obtained  by  equating  the  sample  mean 
20 

X  ■  —  I  X^*'  to  E{X*|0,  t  *  oo }  and  solving  for  0  .  In  other 
m  m  .  n,r  1  o 

words,  if  t 


00  then 
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is  an  unbiased  estimator  of  6  ,  having  a  variance 


/v  7  r-1 

(2.2)  V{0(X  )|6#  t  -  -  0  I 

m  1  o  .  / 

j*o  (n 


In  the  present  section  we  generalize  this  estimator  for  the  case 
of  0  <  tQ  <  00  and  determine  its  asymptotic  properties.  Since  t  is 

fixed,  we  simplify  the  notation  by  letting  E(0)  =  E{X  |0,  tQ} 

*  * 

Since  the  sample  represents  m  i.i.d.  random  variables  X-  , . . .  ,X  , 

1  m 

— *  1  m  * 

the  sample  mean  X  «  —  Z  X.  is  an  unbiased,  strongly  consistent 
■  m  j-i  j 

A  «Jl 

estimator  of  E(6)  .  Accordingly  we  say  that  6  (X  )  is  a  moment 

m 

equation  estimator  (MEE)  of  6  if  it  is  the  root  of  the  equation 


r-1 


(2.3)  6  E  n— T  n,  l-e’S/0)]  =  X*  . 

j*o  J 


The  left-hand  side  of  (2.3)  is  the  function  E(6)  ,  which  is  defined 
on  the  domain  (O,00)  . 

Lemma  2.1 

(i)  E(0)  lj3  strictly  increasing  function  on  (O,00)  ; 

(ii)  E(0)  is  concave  on  (1,°°)  ; 

(iii)  lim  E(0)  =  t  . 

0-xn  ° 

Proof : 

(i)  straightforward  differentiation  yields 
(2.4)  ~  B(j;  n,  p)  *  b^;  n»  P>  * 

and  every  0  <  p  <  1  ,  where  b(j;  n,  p)  is  the  Binomial  p.d.f.,  i.e., 
b (j ;  n,  p)  -  B(j ;  n,  p)  -  B(j-1;  n,  p)  . 

-  8  - 


(2.5)  Jq  E(0) 


I  [1-B(j ;  n,  1-e  °  )] 

.  n— j 

j=o  J 


-t/0 


t  r-l 

*-f  2  b(j;  n,  1-e  °  ) 

j-o 


r-l 


l  -t  /Q 

I  — r  (1-B(j;  n,  1-e  )] 

j»o  n_J 


B(r-1;  n,  1-e 


Comparing  (2.5)  with  (1.9)  -  (1.10)  we  obtain  that 

t 


(2.6)  ^  E(0)  =  px(-|)  >0  .for  all  0  in  (0,®)  . 


This  proves  that  E(0)  is  a  strictly  increasing  function  over  (0,®)  . 


(ii)  Differentiating  (2.5)  we  obtain 


.2  t  -t  /0  t 

»  IT*  /  fl  \  _  O  1*  / _  1  .  _  1  O  V  .  O  fk  / _  «  •  O 


-t/0 


(2.7)  — ~  E(0)  = - x  B(r-1;  n,  1-e  )  +  —x  B(r-1;  n,  1-e 


d0 


e‘ 


(n-r+l)t 


b(r-l;  n,  1-e  °  *  * 


Thus,  for  every  0  >  1  ,  — x  E(0)  <  0  ,  which  proves  the  concavity 

d0Z 


of  E(0)  over  the  interval  (1,®)  . 


r-l  ,  -t/0 

(2.8)  lim  E(0)  -  lim  0  Z  [1-B(j;  n,  1-e  °  )] 

0-H»  0-KO  11  J 


Starting  with  j  *  0  we  obtain 


-t/0 

(2.9)  lim  0  [1-B(0;  n,  1-e  )] 

0-H» 


-nt  /0 


lim  0(l-e 


)  =  nt 


VV 
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On  the  other  hand,  for  every  j  >  1  , 

-t  /e  -nt  /e 

(2.10)  lim  eil-B(j;  n,  1-e  )]  «  lim  0(l-e  °  J 

6-+00  Q-*co 

j  n  _tQ/9  i  ~(n-i)t/6 

-  lim  l  (")  (1-e  °  )i  e  °  -  0  . 

e-H»  i«i  1 


Indeed, 


-t  /e  .  -(n-i)t  /e 

(2.11)  lim  0(l-e  °  )A  e  ° 


,  if  i-1 


,  if  i  >  1 


Thus,  from  (2.8),  (2.9)  and  (2.10)  we  obtain  the  result. 


Corollaries : 


(i)  There  exists  a  unique  solution  to  equation 


(2.3)  ,  0(X*)  ,  for  each  0  <  X*  <  t  . 

m  mo 


(ii)  ®(X*)  does  not  exist  if  X*  =  t  . 
m  mo 


(Q.E.D.) 


The  solution  of  equation  (2.3)  can  follow  the  Newton-Raphson 
iterative  procedure 


(2.i2)  e(i+1)  =  §(1)  -  -E(- _ Li 

E"(6(i)) 


,  i-1,2,... 


where  6  is  any  initial  solution,  and  E'(0)  is  the  derivative  of 
E(0)  ,  given  by  (2.5).  If  Mg  denotes  the  number  of  censored  observations 

in  a  sample  of  m  ,  then  P  *  M_/“  is  a  strongly  consistent  estimator 
-to/0  3  s 

of  B(r-1;  n,  1-e  )  ,  as  .  Thus,  a  solution.  0  .  of 


(2.13)  Po  -  B(r-1;  n,  1-e  ~  ) 

8 


is  a  consistent  estimator  of  6  and  could  serve  as  an  initial  value 
*  / 1 ) 

0in  the  sequence  defined  by  (2.12).  If  M  ■  0  we  can  start  with 

r  1  ^ 

the  "unbiased"  estimator  0  «  X*/  I  ■  \  . 

"  j-o  n-J 

A 

We  derive  now  formulae  for  the  asymptotic  mean  and  variance  of  0(X  )  , 

D 

as  m-*»  . 

According  to  the  Central  Limit  Theorem,  \Zm~(X*  -  E(0))-^ 

m 

2  *1 

N(0,  o  (0;  tQ))  .  bet  E  (x)  denote  the  unique  inverse  of  E(0)  . 

Notice  that  the  MEE  is  the  value  of  E  *(X*)  .  Moreover,  since 

m 

E  1(E{X*  ;  0,  tQ}  )  *  E  *(E(0))  *  0  ,  we  consider  the  Taylor  expansion 


of  E  (x)  around  E(0)  and  obtain 


(2.14)  0(X*)  -  0  +  (X*  -  E(0))  — —  +  <>  (-^~)  ,  as  r». 

m  m  E^(0)  p  /m 


Thus, 


(2.15)  Jm~  (0(X*)  -  0) 
m 


(  A0:  to)  \ 

N  0,  - r 

\  (E '(Q)r  / 


»  a • s • 

Moreover,  since  X  - *  E(0)  ,  and  since  E  (x)  is  a  continuous 

m 

transformation,  0(X*)  is  a  strongly  consistent  estimator  of  0  ,  as 

m 

A  _  a 

m-»«>  .  The  asymptotic  variance  of  0(X  )  is  thus, 

m 


0  (0;  t  ) 

(2.16)  AV{0(X*) ;  0,  t  }  =  - 

"  °  m(E'(8))2 


This  asymptotic  variance  will  be  compared  later  with  the 
asymptotic  variance  of  the  MLE,  to  determine  the  asymptotic  relative 
efficiency  of  the  MEE. 


In  order  to  obtain  an  expression  for  the  asymptotic  bias  of  the 
MEE,  we  add  another  term  to  the  expansion  (2.14). 

The  second  order  derivative  of  E  *(x)  is 


(2.17) 


d_ 

dx 


JL _ 

E'*(E_1(x)) 


E<y(E~1(x)  ) 
(E'(E_1(x)))3 


Hence,  we  obtain  the  expansion 


(2.18)  0(X*)  -  0  +  (X*  -  0) 


m 


E  (0) 


-  \  (X*  -  0) 2  — ^  3  +  0  (i)  , 

2  “  (E"  (0) ) 3  P  m 


A 

nr*00  .  It  follows  that  the  asymptotic  bias  of  0(X  )  is 

m 


(2.19)  B(0;  tQ)  = 


-o  (0;  to) 
2m 


E"(0) 


(E'(0))‘ 


According  to  Lemma  2.1  (ii)  ,  B(0;  t  )  >  0  for  all  0  >  1  . 

o  - 


3.  Maximum  Likelihood  Estimation 

Let  X*  =  (X?  ,...,X*)  be  a  vector  of 
cj&  l  in 

The  likelihood  function  of  0  based  on  X*  , 

r\J& 

proportionality,  is  given  by 

*  m  *  1i 

(3.1)  L(0;  X*)  -  n  [f  ^(X. ;0) ]  [B( 

'um  n»r  i 


m  i.i.d.  random  variables 
up  to  a  factor  of 


t  1-1 
r-1;  n,  F(^))] 


where 

,  if  X*  <  t 
i  o 

,  if  X*  ■  t 
i  o 

■El.  and  let  £,(0;  X*)  denote  the  log-likelihood 

1  ».  Ill 


(3 


-  -■  s: 


Let  K 


function,  then 


t 

(3.3)  f.(0;  X*)  -  -  K  log  6  +  (m-K  )  log  B(r-1;  n,  F(-f)) 

m  o  y 

»  * 

+  Z  I  [(r-1)  log  (l-exp{-  X  /0})  - 

i-1 

(n-r+1)  xJ/0] 

The  maximum  likelihood  estimator  of  0  (MLE)  is  the  value  of  0 
in  (O,00)  which  maximizes  Jl(0;  X*)  . 

Consider  the  score-function,  S(0;  X*)  ,  which  is  the  partial 
derivative  of  £(0;  X*)  ,  with  respect  to  0  .  According  to  (2.4)  we 
obtain  that 


K  t 

(3.4)  S(0;  X*)  =  -  — jr  +  (m-K  )  (n-r+1)  ~  . 

rja  v  m  gZ 

-t /e 

b(r-l;  n.  1-e  ) 

B(r-1;  n,  1-e  to^6) 

m  (n-r+1 )X* 


Lemma  3.1.  If  K  =  0  then  there  is  no  MLE  of_  0  . 

Proof:  If  I  *  0  for  all  i  then,  according  to  (3.4) 
S(0 ;  X*)  >  0  for  all  0  in  (0,“)  . 


-X*/0 

,  X*  e  1 

(r-1)  .  _ i _ 

02  1-exp {-X*/0} 


Theorem  3.1.  If  K_  >  1  then  there  exists  a  unique  MLE 

|  Proof 

• 

• 

The  equation  S(0;  X  )  =  0  is  equivalent  to  the  equation 

r\jn 


(3.5)  6 


_n _  ?  *  _j£l  "  h  Xi 

Km  i«l  1  i_  Km  i-1  =Si75’ 

i-e 


+  ^  <-V  <«♦!> 

"  B(r-li  n,  1-e  °  ) 


Let  H(0;  X*)  denote  the  R.H.S.  of  (3.5).  We  show  first  that,  for 

r\jm 


fixed  X*  , 

rjn 


(i)  lim  H(6 
0-»o 


^  < + *•  H 


(ii)  lim  H(0;  X*  )  -  -  •  . 

e-*»  m 


Indeed, 

(3.6) 


-t  /0 

b(r-l;  n,  1-e  °  ) 

**t  /0 

B(r-1;  n,  1-e  ) 


r-2  (*)  /  t  /0  V(r-l-j)"|-l 

1  +  E  -i-4-  le  °  -l) 

[  J-  U-i)  V  '  J 


Thus, 


(3.7) 


-t  /e 

b(r-l;  n,  1-e  )  _ 

-t  /e  " 

B(r-1;  n,  1-e  ) 

-t  /0 

b(r-l;  n,  1-e  °  ) 

-t  /0 

B(r-1;  n,  1-e  °  ) 


rsrw 
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From  (3.5)  and  (3.7)  we  obtain  (i)  and  (ii).  We  show  now  that 

■go”  H(0;  X^)  <  0  for  all  0  in  (0,°°)  .  This  will  prove  that  there 
is  a  unique  finite  0  satisfying  (3.5). 


Indeed, 


(3.8)  H(0;  j£) 


r-1 


2  ^  *i 

k  ez  i*i  1 


m 


*2  -*V6 

X*  e  1 

~Xi/B  2 
(1-e  1  )Z 


+  r  <»V  ("-r+1>  -je 

IQ 


1  ♦  ?  fil 

(r-l)  1  ' 


But 


(3.9) 


90 


[,.s  a 


t  t  /  0r  r-2 
o  o 

~2  e 
0Z 


-2 


r-2  (?)  /  t  /e  \-(r 

'  £  (A)  <r'1'J)  (e  "  ' 


-(r-j) 


substituting  (3.9)  in  (3.8)  we  obtain  that 
for  all  0  in  (0,<»)  . 


9 

90 


H(0;  X*)  <  0 


(Q.E.D.) 


Numerical  examples  have  illustrated  that  the  function  H(0;  X*) 

might  be  very  steep  in  the  neighborhood  of  0  ,  and  the  Newton-Raphson 
method  has  proven  to  be  instable.  We  have  used  therefore  a  search 
procedure  for  the  maximum  of  the  log-likelihood  function,  which 
provided  good  numerical  results. 
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4.  The  Fisher  Information  Function  and  the  Asymptotic 
Distribution  of  the  MLE 

The  time-censored  r-th  order  statistic  X*  has  a  c.d.f. 

G  (y;  6,  t  )  ,  which  is  specified  in  (1.2).  The  corresponding 
n,r  o 

generalized  p.d.f.  is 

(4.1)  g  (y;  0,  t  ) 

11)1  U 

We  observe  that  these  p.d.f's  have  the  following  properties: 

(i)  The  support  of  g  (u;  0,  t  )  is  (0,t  ]  , 

n,r  o  o 

independently  of  6  . 

(ii)  gR  r(y;  0,  tQ)  has  continuous  partial  derivatives, 

with  respect  to  0  ,  at  every  ye  [0,  tQ]  .  Moreover,  since 
their  support  is  independent  of  0  , 

- 5 —  8_  _  (y;  0»  t  )  is  a  uniformly  integrable  function  of 

36  n,r  ° 

y  ,  for  every  isl,2,...  and  all  0  . 

(iii)  The  p.d.f's  satisfy: 

t  t 

W’2>  “tr{  i  Gn.r<y;  6>  V  ‘  {  if-  tdGn,t(yi  6’  'o’1  *  0 

for  all  0  e  (0,°°)  .  This  is  proven  by  the  following 
Lemma  4.1 

(4.3)  Eeyj-  log  gR  r(X*;  0,  tQ)  J  -  0  ,  all  0e(O,o°)  . 

Proof: 


J  £n.r  <f> 


»  o  <  y  <  t 


B(r-1;  n,  F(-^))  ,  y  «  t 


As  in  (3.3)  we  obtain 


(4*4)  “ST  lo«  «n,r(y;  e*  to) 


1  r-1 


e  e2  i_e-y/e 


n-r+l  _  . 

+  - —  y  ,  0  <  y  <  t 

0^  o 


i  (n-m,  mat  - 

B(r-1;  n,  1-e  0  ) 


*  y  *  ‘o  ‘ 


Thus, 


(4.5)  EQ{-^-  log  g„  ,(X*;  6,  tj) 


|  r  “t  /9  ^  t  - 

-  q-  n,  1-e  °  )J  +  -j  (n-r+l)  b(r-l;  n,  1-e 


n~r-t1  /  y  f  ,(y;0)dy  + 
ez  o  n,r"1 


Furthermore,  according  to  (1.9), 


S  y  <n,r<ywy 

0  o 


(4-6)  Ee  Hr10*  *n,r«*'  6’  '„>l 


i  r  -t_/e  -i  t 

-  £  j^l-B(r-l;  n,  1-e  )J  +  (n-r+l)  b(r-l;  n,  1- 

,  n-r+l  I  /  ,  ,  ,  .  -x/0,.  ,  .  ,  .  .  *to/e 

+ - 5~  I  /  b(r-l;  n,  1-e  )  dx  -  t  b(r-l;  n,  1-e 

0*  L  o  ° 


)  dx  -  t  b(r-l;  n,  1-e 


Finally, 


(4.7)  /  b(r-l;  n,  l-«f*/0)  dx  -  (^)  f°  (l-e'*'8)'"1  . 


.e-(n-r+1)x/e  dx  -  jJjj  1  d, 

a  r  -t  /Q  t 

•  TT-  1-B(r-1;  n,  1-e  °  )  . 


Substituting  (4.7)  into  (4.6)  we  obtain  (4.3). 


(Q.E.D.) 

According  to  Lemma  4.1,  the  Fisher  Information  Function  is 

(4.8)  1(8:  t0)  -  Eo  {(  -jf-  log  go  r(X*i  8.  t0)]2}  . 
Equivalently, 

t  « 

(4.9)  1(0;  to)  -  Q(-f)/Q 


where 


-y 

(4.10)  Q(n)  -  f  £  _(y)  U+(r-l)  ^ 

o  n’r  1-e' 


(n-r+1)  y]2  dy  +  n2 (n-r+1)2-  b  illlx  n».  1-e  ^ 


B(r-1;  n,  l-e”n) 


Expansion  of  the  quadratic  form  in  (4.10)  and  some  algebraic  manipulations 
yield  the  formula 


(4.11)  Q(n)  -  -  [1-B(r-1;  n,  l-e~n)] 

+  2n  (n-r+l)  b(r-l;  n,  l-e”1^) 

+  r)2  (n-r+1) 2  b2(r-l;  n,  l-e”r*)/B(r-l;  n,  l-e”71) 
n-r+1  (2) 


+  n  K‘t-2(T]) 


+  (n-r+1)2  A2  p<2>(n)  . 


(2), 


where  the  second  incomplete  moment  '  (r))  is  given  in  formula  (1.11) 


It  is  interesting  to  study  the  effect  of  the  time-censoring  on 
the  Fisher  Information  Function.  For  this  purpose  we  prove  that  when 


t  ■  00  then 
o 


(4.12)  1(6)  -  I  (6 ;°°)  -  Q/6  , 


where 


>.13)  Q  -  /  £  (y)  [l+(r-l)  -?-■  -  -  (n-r+l)y]2dy 

r  -  -y 

o  l-e 


■  1  +  (r-1)2  f  y2  — Cl 


fn,r«  ** 


+  (n-r+1)2  /  y2  f  (y)  dy 

n*  1 


00  — y 

r  e  J 


+  2(r-l)  /  y  - — -  f  (y)  dy 

1“y  ll  4 1 

-e  * 


-  2  (n-r+1)  /  y  f  (y)  dy 
n  n»r 


-  2(r-l)  (n-r+1)  /  y2  - - —  f  (y)  dy  . 

o  l-e  y  n,r 


Substituting 


(4.14) 


— * - f  (y) 

1-e-y  n*r 


n-r+1  ,  .  . 

r-1  fn,r-l^y^  ’ 


(4.15) 


■e.  —  .  f  (y)  M  1  (n~r~>~2 )  £  /y\ 

(1-e-y}2  n,rW  (r-1)  (r-2)  n,r-2W 
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and  applying  the  formulae 


(4.16) 


/  y2  {n  <(y)  ■  1  — — 2  *  (  1  niiti  )2  > 

o  i-1  (n-i+1)2  Vi-l  n  1+1  / 


we  obtain  after  some  manipulations 

(4.17)  Q  »  1  +  2  i  -1—  +  2 

n-r+2  i,1  n-i+1  (n_r+2)2 


P(n-r+l)  f  r‘2  JL _ +  /  r~2  _1 _ \21 

r"2  Li-1  (n-i+1)2  \i*l  n“i+1  /  J 


As  shown  in  the  present  section,  all  the  conditions  for  the 

consistency  and  asymptotic  normality  of  the  MLE,  §  ,  are  satisfied 

m 

(see  Lehmann  [4]).  Thus,  for  a  fixed  6  , 


(4.18)  0  -  0  -  0  (  ,  as  nr*00 

1,1  P  vm 


(4.19) 


The  convergence,  however,  is  not  uniform.  For  values  of  0  close  to  t( 
the  probability  Pq{K  -  0}  could  be  close  to  1.  For  example,  if 

U  ID 

0  -  tQ  ,  n-30,  r-27  ,  B(26;  30,  1-e"1)  -  .999  . 

Thus,  in  a  sample  of  size  m-50  ,  P{K  -  0}  -  (.999)^  «  .9512  .  Thus, 

ID 

the  probability  is  over  .95  that  the  MLE  would  not  exist.  But,  in  a 
sample  of  size  m  -  1,000  this  probability  drops  to  .368  . 
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In  the  present  section  we  numerically  illustrate  the  asymptotic 
variances  of  the  MLE  and  MEE,  and  provide  some  simulation  results.  These 
simulations  demonstrate  the  actual  behavior  of  these  estimators  in  samples 
of  size  m*20, 30,40  and  50.  The  numerical  computations  of  the  present 
section  are  restricted  to  the  case  of  n=30  and  r*27  .  These  are 
the  parameters  used  by  Epstein  et  al  [2]  in  the  photodynamic  bioassays. 

The  MLE  estimator  6  and  the  MEE  estimator  0(X*)  are  both 

m  m 

scale  equi variant  in  the  sense  that,  if  X  -*■  0  X  for  some  0  >  0  , 

n>  r  n  *  r 

and  correspondingly  t  -*■  0  t  then  0  +00  and  0(0X*)  *  0  §(X*)  . 

o  o  mm  m  m 

Moreover,  since  p  =  t  /0  is  invariant  under  such  transformations,  the 

o 

(asymptotic)  variances  of  0  and  0(X  )  under  0  are  equal  to  0 

m  m 

times  the  corresponding  variances  under  0  *  1  .  We  therefore  report 
here  the  variances  of  the  estimators  under  the  case  of  0  =  1  and 

different  values  of  p  =  tQ/0  •  In  Table  5.1  we  provide  values  of 

Q(p)  for  some  values  of  p  .  This  function  was  determined  according 
to  (4.11)  for  0  <  p  <  00  and  according  to  (4.17)  for  p  =  <» 

(uncensored) .  Notice  that  Q(p)  is  the  value  of  the  Fisher  information 
function  at  0  *  1  ,  i.e.,  1(1; p)  -  Q(p)  . 

Table  5.1 

Values  of  Q(0)  for  n*30,  r»27 
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As  seen  in  Table  5.1,  the  amount  of  information  per  observation 
is  maximal  when  there  is  no  censoring  (p*00)  .  Moreover,  if  p  >  3 

the  amount  of  information  is  close  to  the  maximal  one.  The  asymptotic 
efficiency  of  the  MLE,  under  p  ,  compared  to  the  no-censoring  case  is 

(5.1)  AE(p)  -  Q(p)/Q 

Thus,  as  seen  in  Table  5.1,  AE(1.8)  *  .564  .  On  the  other  hand, 

AE(3)  *  .99  .  Thus,  for  6  values  close  to  the  censoring  point  t 

we  need  considerably  larger  samples,  than  in  the  uncensored  case,  to 
attain  a  specified  asymptotic  precision. 

In  Table  5.2  we  present  the  values  of  the  asymptotic  variances 
of  the  MEE  and  the  MLE,  which  have  been  determined  according  to  formulae 
(2.16)  and  (4.19),  respectively.  The  asymptotic  relative  efficiency 
(ARE) 

Table  5.2 

The  Asymptotic  Variances  of  the  MEE 

_3 

and  the  MLE  [in  units  of  10  ],  for  0*1,  m=50 

n=30,  r=27  and  the  ARE 


a 

MEE 

MLE 

ARE 

1.80 

2.2723 

1.8819 

.828 

2.00 

1.6549 

1.4574 

.881 

2.25 

1.3153 

1.2269 

.933 

2.57 

1.1450 

1.1152 

.974 

3.00 

1.0809 

1.0729 

.993 

of  the  MEE  relative  to  the  MLE  is  defined  as  the  ratio  of  the  asymptotic 
variances,  i.e., 

AV{0  ;  1,  p} 

(5.2)  ARE(p)  -  - r2— - 

AV{0 (X*) ;  1,  P> 

ID 
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As  seen  in  Table  5.2,  as  p  increases  the  value  of  the  ARE  approaches 

one.  In  other  words,  the  asymptotic  efficiency  of  the  MEE  is  close  to 

that  of  the  MLE  when  the  value  of  6  is  considerably  smaller  than  that 

of  t  ,  i.e.,  p  >  3  .  If  we  suspect  that  6  is  close  to  t  the 
o  -  o 

MLE  is  a  more  precise  estimator  than  the  MEE. 

It  is  also  interesting  to  study  how  the  two  estimators  perform 

when  the  samples  are  not  very  large.  For  this  purpose  we  have  conducted 

a  series  of  simulations.  In  these  simulations  we  generated  random 

* 

samples  of  m  values  of  X  ,  for  n  =  30  ,  r  =  27  ,  t  =90  and 

n,r  o 

0  =  30  ,  35,  40,  45  and  50  .  For  the  sample  size,  m  ,  we  considered 
the  values  m  «  20  ,  30,  40  and  50  .  For  each  of  these  samples  we 

A  —  JL  A 

computed  the  values  of  the  MEE,  9(X  )  ,  and  of  the  MLE,  6  .  These 

m  m 

trials  were  repeated  M=500  times  independently.  If  all  the  values  of 

X*  were  censored,  a  message  "NO-EST."  was  printed.  It  is  interesting 
n,  r 

to  notice  that  in  the  case  of  m=20  and  6=50  (p=1.8)  we  obtained  this 

message  two  samples  out  of  500  independent  samples.  This  conforms  to  the 

theory,  since  in  the  case  of  p=1.8  ,  the  probability  of  censoring  is 

B(26;  30,  .8347)  =  .7543  .  Thus,  the  probability  that  all  m=20 

20 

observations  will  be  censored  is  (.7543)  =  .0036  .  If  such  samples 

are  simulated  independently  M=500  times,  the  expected  number  of  samples 
with  complete  censoring  is  1.8  . 

In  Table  5.3  we  present  the  means  and  mean  square-errors  of  the 
simulated  values  of  the  MEE  and  MLE.  It  is  interesting  that  in  all  cases 
the  bias  is  negligible.  In  addition,  the  MSE's  of  the  MEE  are  not 
significantly  larger  than  those  of  the  MLE  when  6  <  40  (p  >  2.25)  . 
Furthermore,  when  m=50  the  MSE's  of  the  MLE  and  of  the  MEE  are  not 
significantly  larger  than  the  corresponding  asymptotic  variances,  even 
if  6  >  40  (p  <  2.25)  .  For  example,  when  m=50,  0=50  ,  p=1.8  and 
from  Table  5.2,  the  asymptotic  variance  of  the  MLE  is  2.5  x  1.8819  =  4.705. 
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Table  5.3 


Means  and  MSE  of  M-500  simulated 


values  of  the  MEE  and  MLE,  for  n-30,  r«27 


t  -90,  0-30(5)50,  and  m-20(10)50 


30.25 

35.25 
40.36 
45.49 
50.54 


2.358 

3.445 

5.475 

8.444 

18.938 


30.21 

35.21 
40.28 
45.37 
50.50 


1.779 

2.520 

3.854 

5.824 

9.400 


30.27 

35.28 
40.30 
45.33 
50.41 


1.325 

1.866 

2.642 

3.943 

6.972 


30.29 
35.31 

40.30 
45.35 
50.43 


-The  estimates  are  based  on  M  -498  samples 
Two  samples  were  completely  censored. 
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29.97 

35.02 

40.11 

45.30 

50.57 


29.94 

34.96 

40.03 

45.16 

50.47 


29.98 

35.00 

40.03 

45.10 

50.33 


30.02 

35.02 

40.04 

45.10 

50.27 


2.319 

3.606 

5.789 

9.795 

23.680 


1.71 

2.57 

4.08 

6.53 


12.944 


1.223 

1.793 

2.693 

4.470 

8.586 


0.948 

1.364 

2.082 

3.382 

6.023 


The  corresponding  standard-error  of  the  variance  estimate  is  approximately 

SE*4.705  (2/500) 1/2  -  .298  .  Thus,  4.705  +  2  SE  -  5.300  ,  which  is 
above  the  MSE  value  of  the  MLE  in  Table  5.3  .  Similarly,  the  MSE  of 
the  MEE,  for  m*50,  0*50  ,  is  within  the  . 95-confidence  interval 
constructed  around  the  corresponding  asymptotic  variance.  This  is  not, 
however,  the  case,  when  m*20  and  0  >  45  .  The  MSB’s  of  the  MLE  and  the 
MEE  are  significantly  larger  than  the  corresponding  asymptotic  variances. 
For  m*20  ,  0*50  ,  the  asymptotic  variance  of  the  MLE  is  11.762  .  On  the 
other  hand,  the  estimated  MSE  is  18.938  .  This  shows  that  the  actual 
variances  of  the  MLE  and  the  MEE  are  considerably  larger  than  the 
asymptotic  variances  if  the  samples  are  not  sufficiently  large  and 

p  <  2.00  . 
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